Light Higgs Mass Bound in SUSY Left-Right Models 
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We show that in supersymmetric left-right models (SUSYLR), the upper bound on the lightest neutral Higgs 
mass can be appreciably higher than that in minimal supersymmetric standard model (MSSM). The exact mag- 
nitude of the bound depends on the scale of parity restoration and can be 10-20 GeV above the MSSM bound if 
mass of the right-handed Wr is in the TeV range. An important implication of our result is that since SUSYLR 
models provide a simple realization of seesaw mechanism for neutrino masses, measurement of the Higgs boson 
mass could provide an independent probe of a low seesaw scale. 



1. Introduction One of the main missing links of the 
otherwise immensely successful Standard Model (SM) is the 
Higgs boson which plays the crucial role in giving masses to 
all elementary particles in nature. It is therefore rightly the 
focus of a great deal of theoretical |1] and experimental en- 
quiries. Even though the Higgs boson mass in SM is arbi- 
trary, some ideas about how heavy the Higgs boson can be 
gained in the context of different plausible extensions of SM 
as well as from other considerations [2^3]. Typical upper lim- 
its from, say unitarity considerations 1131 is in the TeV range. 
This bound is however considerably strengthened in one of 
the most widely discussed possibility for TeV scale physics, 
supersymmetry. Specifically in the minimal supersymmet- 
ric SM (MSSM), the upper bound on the Higgs boson mass 



is M™°^ < 135 GeV ||4|] when one and two loop radiative 
coiTections are included. Present collider searches provide a 
lower bound on the SM Higgs mass L5J of 114 GeV leaving a 
narrow region which need to be probed to test MSSM. If the 
Higgs mass is found to be above this upper limit, does it mean 
that supersymmetry is not relevant for physics at TeV scale? 
The answer is of course "No" since there exist simple and 
well motivated extensions of MSSM, e.g. the next-to-MSSM, 
which extends the MSSM only by the addition of a singlet 
field la] where there is a relaxation of this bound to about 142 
GeV or so |7[l- There are also other examples in literature|8] 
where simple modifications of the post-MSSM physics can 
provide additional room for Higgs mass. 

In this paper we discuss an alternative scenario motivated 
by neutrino mass as well as understanding of the origin of par- 
ity violation |9] where the upper limit on the light Higgs mass 
is relaxed compared to MSSM. The model is the supersym- 
metric left-right model (SUSYLR) ^ [H with TeV scale 
parity restoration (or TeV right-handed gauge boson mass 
Wr)- The change in the Higgs mass upper limit comes from 
the contribution of the D-terms and satisfies the decoupling 
theorem i.e. as the Wr mass goes to infinity, the Higgs mass 
upper bound coincides with that for MSSM. This effect is to 
be expected on general grounds ||13[1 in gauge extensions of 
MSSM. 

2. Basics of the SUSYLR model The gauge group of 

this model is SU{2)l x SU{2)r x U{1)b-l x SU{i)c. The 
chiral left-handed and right-handed quark superfields are de- 
noted by Q = {u,d) and Q"^ = [u'^^d'^) respectively and 



similarly the lepton superfields are given by L = {v, e) and 
IjC ^ ^-^c^ gc-j 'pjjg Q ^jjj ^ transform as left-handed doublets 
with the obvious values for the B — L and the Q"^ and L'^ trans- 
form as the right-handed doublets with opposite B — L val- 
ues. The symmetry breaking is achieved by the following set 
o_f Higgs superfields: $„(2, 2, 0, 1) (a = 1, 2); A(3, 1, +2, 1); 
A(3, 1, -2, 1); A'=(l, 3, -2, 1) and A=(l, 3, +2, 1). We in- 
clude a gauge singlet superfield S to facilitate the right handed 
symmetry breaking. The symmetry breaking can also be car- 
ried out hy B — L = 1 doublet fields, for which our results 
also apply. A virtue of using triplet Higgs fields is that they 
lead to the see-saw mechanism for small neutrino masses us- 
ing only renormalizable couplings. In addition, as was noted 
many years ago [11], low scale Wr^ requires that R-parity 
must break spontaneously. This leads to many interesting phe- 
nomenological implications that we do not address here. 
The superpotential for the model is given by: 



W 



S* [Tr (AA + A^A'^) - 4] (1) 



In order to analyze the Higgs mass spectrum, we write down 
the Higgs potential for the model including the soft SUSY- 
breaking terms: 
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where Vp and Vd are the standard F-term and D-term po- 
tential and Vs is the soft-SUSY-breaking terms which can be 
found in the literature 111. 121 . Minimization of the Higgs 
potential leads to the following vacuum configuration for the 
A'^ and ly" Higgs fieldslflTll: 
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Note that in the SUSY limit vr = vr and (i^=) = 0. _ 
the presence of supersymmetry breaking terms however, (z/'=) 
is nonzero. On the other hand, if this model is extended to 
include a B-L=0 right handed triplet with nonzero vev, there 
appears a global minimum of the potential which has {f'^) = 
UM even in the presence of susy breaking terms. Since the 
vev of ly'^ is not relevant to our discussion, we will work with 



B-L=0 triplet model and set (1/^=) = henceforth. The SM 
symmetry remains unbroken at this stage and is broken by the 
vevs of the $ fields. We can write these fields in terms of their 
MSSM Higgs content: 
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with vevs {H', 

Before proceeding to discuss upper bound on the light neu- 
tral Higgs mass in this model, we wish to make a few com- 
ments on the implications of the TeV scale Wr models for 
neutrino masses. First, in the non-SUSY left-right model 
where neutrino mass has both type I and type II seesaw con- 
tributions, having a TeV scale Wb is unnatural since the 
type II seesaw contribution then becomes extremely large. 
This is due to the presence of non-zero couplings of type 
TiipALfj)^ A'^jj, which are allowed by the symmetries of the 
theory. On the other hand, in the SUSYLR model, this cou- 
pling is absent due to supersymmetry and therefore there is 
no type II seesaw contribution to neutrino mass. As far as the 
type I contribution is concerned, if we choose h^ ~ h^ where 
hi, is the Dirac neutrino Yukawa coupling, then we can have a 
few TeV Wr and neutrino masses of order of eV Thus as far 
as neutrino masses go, low scale Wr is a realistic model. 



M^± — Qrv]^ because the triplet vev breaks custodial sym- 
metry for the right-handed sector. 

Since there is no coupling between the bidoublet Higgs and 
the triplet Higgs fields responsible for parity breaking in Eq. 
(1), any change in the effective MSSM doublet Higgs poten- 
tial below the vr scale must originate from the D-terms, 

Vd = 4|Tr[2A=V™A= + 2A^t^^A^ + ci>r^ci>t]P 



^ (TrpA'^^A'^ - 2A'=^A'=])^ 



The contribution to the neutral Higgs fields couplings is. 



(4) 



ynent. ^ ^ I Xr[$r3^$t] I 



g%VRRe[{A^° - A=")]Tr[$r3^<l>t] 



^Ij + sIlH 



A 



cO 



V2 



(5) 



The coupling is linear in the field Re[(A'^° — A''")], which 
we will call a^ . As (T_ field becomes heavy, its coupling to 
[$rj$t] will generate new quartic term in the MSSM doublet 
field potential, which in turn will lead to new contributions to 
Higgs mass upper bound. Collecting this new effect, we get 
for the Higgs quartic term: 



3. Light Higgs mass bound: single bi-doublet case We 

proceed to consider the bound on the light neutral Higgs mass 
in the SUSYLR model. We work in the limit where vr and vr 
are much bigger than the SM scale. In this limit, we search for 
additional contributions to the MSSM Higgs potential, which 
will be at the heart of the change in the upper limit of the 
Higgs boson mass. 

We first illustrate this in a one bi-doublet model. This sim- 
ple model leads to vanishing CKM angles at the tree level, 
which can be fixed in one of two ways: (i) by including ra- 
diative correction effects from squark mixings [15] or (ii) by 
including a second bi-doublet which decouples from the low 
energy sector but it has a tadpole induced vev that generates 
the correct CKM angles. We discuss the case (ii) toward the 
end of the paper. The interesting point is that neither of these 
affects the Higgs mass upper bound that we derive. For the 
model under consideration, we first show that in the SUSY 
limit the low energy Higgs potential recovers that of MSSM 
in the same limit. When soft SUSY breaking terms is taken 
into account, there appear new contributions to the MSSM 
Higgs potential, serving to raise the upper limit on the light 
Higgs mass, which can be significant for TeV scale Wr. 

We start with a review of the well known symmetry break- 
ing of the model by the triplet Higgs fields A'^ and A'^ in the 
SUSY case. The gauge bosons get mass from the kinetic terms 
of triplets and after symmetry breaking, the massless gauge 
boson and gaugino corresponding to U{1)y is the combina- 
tion B = f ^ ^ WsR + I /^ ^ Vbl, with the hyper- 

charge gauge coupling given by ^5- = -5- H — 5 — . The 

9y 9r 9bl 

heavy Z'-boson has mass squared M|, = '^{Qr + Ssl)^!;- 
There is a factor 2 compared with the charged Wr boson mass 
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To evaluate this new contribution, we need to know Mcr_. 
This has two potential contributions: (i) from the D-term and 
(ii) from the F-term contribution to the Higgs potential. It 
turns out that in the SUSY limit, the only contribution to Ma_ 
is from the D-terms and we have M^_ = 2{g\ + gBL)''^\- 
This follows not only from actual calculations but also from 
the fact that a^ is a member of the Goldstone supermultiplet, 
all members of which must have the same mass as Z' in the 
SUSY limit. This result would hold even if the superpotential 
had a term of the form jjlA'^A^. Using this in Eq. (|6]l, it is easy 
to see that the net contribution to the quartic term in the Higgs 

superpotential becomes -^ (\Hu\'^ — \Hd\'^) ■ This is noth- 

8 
ing but the Dy contribution to quartic Higgs doublet term in 

MSSM. Since in the decoupling limit, we get MSSM, as ex- 
pected from the decoupling theorem. 

Let us now switch on the supersymmetry breaking terms. In 
their presence, the a^ field has aditional contributions which 
lead to a shift in the Higgs masses. To see this we intro- 
duce soft mass term Trigs'^ S as well as SUSY breaking mass 
terms for the A^ and A'^. Taking the same SUSY breaking 
terms for all the fields gives different value for the cr_ field 
mass and we get for the contribution to the quartic Higgs term 

%^(|i/„p-|i/,P)' where 
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where niQ in the above equation is a generic soft mass term 
for sparticles that breaks supersymmetry. This leads to an 
enhancement of the Higgs mass upper bound since gy efF > 
9y sm ■ To get an idea about how large the change in the 
upper bound is likely to be, we take toq — 1 TeV, vr — 2 
TeV, then from gl ss 0.42 and gy ~ 0.13, we get the ratio 

2 I 2 

r = ^., ^;f^ Ki 1.1 which will give 10% increase of the 
tree level upper bound on the light Higgs mass i.e. it increases 
from Mz = 90 GeV to 100 GeV. 

It is also worth pointing out that as the scale of parity viola- 
tion goes to infinity, this new contribution goes to zero and one 
recovers the MSSM result. This is an important consistency 
check on our result iQi 



4. One-loop radiative corrections and numerical result 

In the following, we discuss the radiative corrections to the 
Higgs boson mass for the model above. It is well known that 
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where the supertrace Str means Str f{M.^) = 
^■(— l)^'''(2Ji + l)f{mf), where the mass m^ is cal- 
culated in background fields, and the sum counts all the 
fermionic and bosonic degrees of freedom. Ji is the spin 
for particle i and Q is the renormalization scale on the order 
of electroweak symmetry breaking. In MSSM, top and stop 
contributes dominantly to AVi because of large Yukawa 
coupling yt, while the bottom and sbottom contribution is 
only important for tan /? ^ 20. 

In SUSYLR model with one bidoublet, we have tan/3 — 
mt/mi, « 40, so the sbottom quark can couples to Hu with a 
large coupling yt. This is from the F-term of bidoublet field 
$1. 

(F*Jy = /jl(QlT2),(r2Q^),+Mll(T2*fT2)j, (9) 

|^$iP = ^inmtcotf] ■VjJr + fiiinifb''j^bR^ 

Note the second term in the second line differs from MSSM 
by a factor ■mt/{nrn, tan/3), since H^ and Hd are unified into 
the same bidoublet. This means the sbottom mass receives a 
large LR mixing proportional to the top Yukawa couplng, even 
for low tan/3 (in the presence of a second bidoublet in the 
realistic model below). Therefore now there are three fields 
that have to be taken into account: top, stop and sbottom. If 
we neglect the small couplings except for yt, and also neglect 
the A-terms their masses can be approximated by 



(10) 



The • • • represents dependence on \Hd\, but since the light- 
est Higgs boson is mainly made up of Hu for large tan/3 
with only a small Hd component, this dependence can be 
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neglected. We can choose a proper scale Q so that the first 
derivative vanishes, which will be important in eliminating the 
explicit Q dependence of Higgs mass, i.e. it can only depend 
on Q through depending on other parameters. Second deriva- 
tive gives radiative corrections to the lightest Higgs mass 
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The first term is as the usual MSSM one. Now we have two 
new terms proportional to /in. Their sum is an even function 
of /ill since changing the sign also interchanges bi ^-> 62- We 
find the sum of second and third term are negative definite for 
arbitrary /in. Actually, for ttiq ^ 1 TeV, one can expand with 
^'-'-™' . The net contribution is non-vanishing only up to the 

2 4 
third order, which is — Shr -tt^^ 



32.^ M'^^mr depending on /in very 
mildly. For /in ranging from 100 GeV to 300 GeV around 
EW scale, this negative contribution is smaller than 1 GeV. 
Note that in this case, we did not have to discuss the details of 
EWSB since it is very similar- to MSSM. 

Let us present the numerical results for the Higgs mass up- 
per bound for this scenario. In Fig. 1, we plot the differ- 
ence in the prediction of upper bound on the lightest Higgs 
boson mass between SUSYLR model and MSSM: A/\/,i = 
^SUSYLR _ j^MSSM p^j. jj^g right-handed scale near 2-3 

TeV, the upward shift of the Higgs mass bound can be of a 
few GeV, increasing as vr decreases. (For symmetry break- 
ing using Higgs doublets, M^„ = {g\ + g%L)v\ and gy.cff 

in Eq. Q gets increased to 4'^j^>"^^l'"l/"l' . Then AMu can 

9R+9BL+"^al "r 

further increase by a factor of 2.) The bound also increases 
with the soft mass scale mo as expected. From the discus- 
sions below Eq. ( fTTI ). non-zero /in always gives small and 
negative contribution. In order not to violate the lower bound 
on chargino mass at LEP2, we choose /in > 100 GeV. 



5. Extension to two bi-doublet case We can extend our 
discussion to the more realistic two bi-doublet case as needed 
to generate the correct CKM angles at the tree level. There are 
two possible ways to do that; both these we discuss below. 
Model A: 

In this case, we identify $1 as the bi-doublet of the previous 
section. We can diagonalize the corresponding Yukawa cou- 
pling matrix h]-,. We are then forced to have all elements of 

the second bi-doublet $2 Yukawa coupling, Kq 7^ 0. Once the 
second $2 has vev, by appropriate choice of this matrix, we 
can generate the desired quark masses and mixings. We will 
see that even though there are four real neutral Higgs fields in 
this case, the upper limit on the light Higgs field remains the 
same as in the single bi-doublet case. 

To see this, first we choose a basis in the ^ space such that 
the superpotential for $'s has the form 



We: 
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FIG. 1 : Higgs mass bound from SUS YLR model shown as the differ- 
ence from MSSM. Tree level results are shown as dashed curves and 
radiative corrected ones are solid curves. Upper, middle and lower 
two curves correspond to mo = 1 TeV (blue), mo = 600 GeV (red), 
mo = 400 GeV (green), respectively. For radiative correction, we 
choose /111 = 100 GeV. For all the curves, we choose tan/3 ~ 40 
as noted. 



We assume that /i22 3> /in. We have then no freedom to 
diagonahze the soft SUSY breaking terms. The sum of the 
Vf + Vs can then in general be written as 



Vf + Vs{<i>^) 



= m^,Tr0|0j + 6iiTr( 
h.c. 



&22Tr( 



l) 



^2 92) 



(13) 



Note that if 777,22 ^ "^12' 'hen the mixed term in the 0's 
will induce a vev for the (f)2 field which is small compared 
to that for the </>! i.e. K2, k,^ ^ 'ti, n'l- To generate the cor- 
rect mass and mixing pattern for the quarks, it is sufficient 
to have the 02 vevs of order of a 100 MeV. For instance if 
mi2 < 10 GeV and 77122 = 1 TeV, then we can estimate 
K2 = ^[11112/17122 = k']^/100 ~ 100 MeV, which is enough 
to generate the strange quark mass as well as other CKM an- 
gles. Note also that one should include the one loop effects 
coming from squark masses and mixings|i5[. While we do 
not give a detailed fit here, it seems clear that this is a realistic 
model where the new Higgs mixing parameter 777,12 is in the 
1-10 GeV range. When it is close to one GeV, the effect on 
the Higgs mass upper bound is also about a GeV lower due to 
off diagonal contributions. We can also keep the $2 Yukawa 
couplings sufficiently small so that their radiative corrections 
do not affect the one loop result. This vacuum then is a per- 
turbation around the vacuum of the single bi-doublet case and 
furthermore due to large /i22, the Hu.d coming from the sec- 
ond bi-doublet will acquire heavy mass and decouple without 
affecting the light Higgs mass upper bound except perhaps a 
small one GeV or so shift. This case corresponds to large 
tan/3 w40. 
Model B: 

In this case, we choose two bi-doublets with the vev pattern 
given by: 
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Yukawa coupling unlike the model A, we can have the value 
of tan/3 = — much lower than 40 by appropriate choice of 
the second Yukawa coupling matrix. There are generally four 
electreweak scale Higgs doublets. Using the standard formula 
in Eq. ([8]l, we have calculated the 1-loop radiative corrections 
to the 4x4 neutral Higgs mass matrix in the presence of SUS Y 
breaking thresholds. As before, we negleted all other cou- 
plings but yt when calculating AVi. We also keep the effect 
of non-zero vev for the right-handed sneutrino. However due 
to small neutrino Dirac Yukawa couplings, the mixing effect 
between the Higgs field and the left-handed sneutrino caused 
by the right-handed sneutrino vev is very small and does not 
affect our result. In order to estimate the upper bound, we 
have done a numerical study to obtain the Higgs mass for ran- 
dom choice of parameters. The results are the scatter points 
in Fig. 2 below for a choice of the generic soft mass scale 
7770 = 1 TeV and right-handed scale vr — 1.5 TeV. Each 
point in the scatter plot represents the lightest Higgs mass for 
a specific choice of parameters. The upper limit therefore cor- 
responds to the topmost set of points in Fig. 2. In contrast, 
the MSSM Higgs mass upper bound is plotted as the yellow 
(lower) curve, which is at most 130 GeV after 1-loop radiative 
corrections with the same choice of ttiq. The red (upper) curve 
is for Model A. We find Fig. 2 that in general SUSYLR model 
the upper bound can be as high as 140 GeV or even more es- 
pecially in the regime 5 < tan /3 < 10. This is higher than the 
prediction of MSSM. Clearly, as the right-handed scale goes 
down, the upper bound inceases. 




Since the down quark masses in this case come from a second 
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FIG. 2: The scatter plot represents the Higgs mass values as a func- 
tion of tan/3 in the two bi-doublet SUSYLR model for random choice 
of parameters. The yellow (lower) curve is the prediction for the light 
Higgs mass upper bound in MSSM, while the red (upper) curve is for 
Model A in SUSYLR. The blue points are for Model B in the SU- 
SYLR case. When plotting the figure, we choose mo = 1 TeV and 
vr — 1.5 TeV. We have used a Monte Carlo simulation for the pa- 
rameter space to generate the plot for Model B. 



6. Comments and Conclusion: Before concluding, we 
wish to make a few comments on the model: 

(i) Low scale non-SUSY left-right models have strong con- 
straints coming from the tree level Higgs contribution to fla- 
vor changing processes. In the SUSY version however, there 
are additional contributions to the same from squark and slep- 
ton sector which can be used to cancel this effect 1, 17,1 . While 



strictly this is not natural, from a phenomenological point of 
view, this makes the model consistent when both the Wr and 
Higgs masses are in the few TeV range. 

(ii) The second point is that unlike other models such as 
NMSSM where the light Higgs mass bound is changed by 
making additional assumptions about the Higgs couplings 
(e.g. not hitting the Landau pole at the GUT scale), in our 
model the increase in the bound is purely gauge coupling in- 
duced and is independent of the Higgs couplings. 

(iii) It is also worth stressing the obvious point that observa- 
tion of a Higgs with mass above the MSSM bound of 1 35 GeV 
is not necessarily an evidence for the SUSYLR model since 
there exist other models with which also relax this bound. One 
needs other direct evidences such as the mass of Wr or Z' 
produced at LHC which when combined with observed higher 
Higgs mass could provide evidence for low left-right seesaw 



scale. 

To conclude, we have pointed out that the upper bound on 
the light Higgs mass is higher if MSSM is assumed to be an 
effective low energy theory of a TeV scale SUSYLR model. 
The increase can be as much as 10 GeV or more depending 
on the scale of parity breaking. If the Higgs boson mass in the 
collider searches is found to exceed the MSSM upper limit of 
135 GeV, one interpretation of that could be in terms of a TeV 
scale seesaw in the context of a SUSYLR model. 

We thank K. Agashe, R Batra, and Z. Chacko for com- 
ments. This work was partially supported by the U. S. De- 
partment of Energy via grant DE-FG02-93ER-40762. R. N. 
M. is supported by NSF grant No. PHY-0652363. Y Z. ac- 
knowledges the hospitality and support from the TQHN group 
at University of Maryland and a partial support from NSFC 
grants 10421503 and 10625521. 



[1] For a review, see S. Dawson, J. Gunion, H. Haber and G. Kane, 
Higgs Hunter's guide, SCIPP-89/13, UCD-89-4, BNL-41644, 
Jun 1989. 404pp. 

[2] for an upper limit from unitarity considerations, see B. W. Lee, 
C. Quigg and H. B. Thacker, Phys. Rev. D 16, 1519 (1977). 

[3] M. Sher, Phys. Kept. 179, 273 (1989). 

[4] H. Haber and R. Hemfling, Phys. Rev. Lett. 66, 1815 (1991); 
J. R. Ellis, G. Ridom and F. Zwirner, Phys. Lett. B 257, 83 
(1991). 

[5] R. Barate et al. Phys. Lett. 565, 61 (2003). 

[6] H. P Nilles, M. Srednicki and D. Wyler, Phys. Lett. B 120, 346 
(1983). 

[7] U. Ellwanger and C. Hugonie, Mod. Phys. Lett. A 22, 1581 
(2007). 

[8] Some examples of such models with relaxed Higgs mass limit 
are: M. Drees, Int. J. Mod. Phys. A 4, 3635 (1989); K. S. Babu, 
I. Gogoladze and C. Kolda, arXiv:hep-ph/0410085; G. Bhat- 
tacharyya, S. K. Majee and A. Raychaudhuri, Nucl. Phys. B 
793, 1 14 (2008); S. Hesselbach, D. J. Miller, G. Moortgat-Pick, 
R. Nevzorov and M. Trusov, Phys. Lett. B 662, 199 (2008); 
M. Dine, N. Seiberg and S. Thomas, Phys. Rev. D 76, 095004 
(2007); V. Barger, P. Langacker, H. S. Lee and G. Shaughnessy, 



Phys. Rev. D 73, 115010 (2006). 
[9] R. N. Mohapatra and J. C. Pad, Phys. Rev. D 11, 566, 2558 

(1975); G. Senjanovic and R. N. Mohapatra, Phys. Rev. D 12, 

1502 (1975). 
[10] M. Cvetic and J. C. Pati, Phys. Lett. B 135, 57 (1984); 

R. M. Francis, M. Frank and C. S. Kalman, Phys. Rev. D 43, 

2369 (1991). 
[11] R. Kuchimanchi and R. N. Mohapatra, Phys. Rev. D 48, 4352 

(1993). 
[12] K. Huitu, P. N. Pandita and K. Puolamaki, arXiv:hep- 

ph/9904388. 
[13] P Batra, A. Delgado, D. E. Kaplan and T. M. P Tail, IHEP 

0402, 043 (2004). 
[14] R. Kuchimanchi and R. N. Mohapatra, Phys. Rev. Lett. 75, 3989 

(1995) 
[15] K. S. Babu, B. Dutta and R. N. Mohapatra, Phys. Rev. D 60, 

095004 (1999). 
[16] The upper bound on light Higgs mass in SUSYLR models de- 
rived in K. Huitu, P. N. Pandita and K. Puolamaki, Phys. Lett. 

B 423, 97 (1998) does not satisfy this decoupling consistency. 
[17] Y. Zhang, H. An and X. d. Ji, arXiv:0710.1454 [hep-ph]. 



